Dynamics of coherent and incoherent emission from an artificial atom in a ID space 
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We study dynamics of an artificial two-level atom in an open ID space by measuring evolution 
of its coherent and incoherent emission. States of the atom - a superconducting flux qubit coupled 
to a transmission line - are fully controlled by resonant excitation microwave pulses. The coherent 
emission - a direct measure of superposition in the atom - exhibits decaying oscillations shifted 
by 7r/2 from oscillations of the incoherent emission, which, in turn, is proportional to the atomic 
population. The emission dynamics provides information about states and properties of the atom. 
By measuring the coherent dynamics, we derive two-time correlation function of fluctuations and, 
using quantum regression formula, reconstruct the incoherent spectrum of the resonance fluorescence 
triplet, which is in a good agreement with the directly measured one. 

PACS numbers: 42.50.Gy, 85.25.-j 



Superconducting circuits with small Josephson junc- 
tions are now commonly used to demonstrate a variety 
of quantum phenomena. The circuits were shown to pos- 
sess discrete energy levels, controlled quantum coherent 
evolution and coherent interaction with quantized 

electromagnetic modes of resonators In a series 

of experiments many fundamental effects from quantum 
optics have been demonstrated fMl^ . Strong interaction 
of the circuits with the electromagnetic fields allows also 
to demonstrate new phenomena, e.g. controllable manip- 
ulations with single photons and photon fields 3, 15| . 

Recently we have realized an analog of a natural atom 
in the open space - an artificial atom in the ID space 
[m . The atom is strongly coupled to a transmission line 
via dipole interaction. Here we study dynamics of coher- 
ent and incoherent dipole emission from such an atom. 
Manipulating the atomic states by microwave pulses, one 
can perform the quantum state and process tomography. 
Differently from the qubit tomography realized by pro- 
jective measurement to its states [l6|, we perform essen- 
tially optical measurements of the dipole emission from 
the non-isolated atom. The excited atom is a dynamical 
system, which continuously emits radiation to the line, 
similarly to natural atoms in the open space. However 
differently from the optical measurements of the natural 
atoms, the collection efficiency of the emission in the ID 
transmission line is very high due to its strong coupling 
to the line llj. We suggest a convenient experimental 
procedure for extracting the two-time correlation func- 
tion of fluctuations, using the time-dependences of the 
coherent dipole emission. The incoherent spectrum of 
resonance fluorescence is then calculated, which shows a 
good agreement with the experimentally measured one. 

As an artificial atom, we use a superconducting flux 
qubit d, inductively coupled to a coplanar trans- 
mission line ^T]\. The two lowest eigenstates of the 
atom (qubit states) are formed by the states associated 
with clockwise and counterclockwise circulating currents 
with the persistent current Ip = 213 nA. Their energies 
are controlled by the external magnetic flux threading 
through the loop $ = $o/2 + (5$, where 5$ is the devia- 



tion from the half-flux quantum $0/2. We tune the flux 
bias to the optimal point (5$ = 0, where the atomic tran- 
sition frequency is uJa/^n = 9.888 GHz. The loop shares 
a part with the transmission line, which results in a mu- 
tual inductance M = 13.6 pH mainly due to the kinetic 
inductance of the shared segment. 

The two-level atom driven by the resonant microwave, 
with the current field Iq cos (ujat — 1^9), is described by the 

CT-e-^'^)/2 in the rotat- 



{ffx ± *o-j^)/2, and 



Hamitonian H = —Ml{a^e'^f 
ing wave approximation. Here cr^ 
ax, CTy, Oz are the Pauli matrices, 51 = 0p/o/(2?i) is the 
Rabi frequency, and (/)p = Ml^ is the atomic dipole tran- 
sition matrix element. The dynamics of such a, system is 
described by a spin-1/2 in the magnetic field |18| and is 
governed by the optical Bloch equations 



where 



B 



— = Ba + 5, 
at 



— r2 —ilsmip 
— r2 —51 COS 93 
flsimp flcosif — Fi 
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(2) 



b — {0,0,— Fi}, and a ~ {{ax), {<Jy), (az)} is a vector 
of the expectation values of the Pauli matrices. This 
vector represents the atomic state, accounting incoher- 
ent processes of relaxation with the rate Fi (decay of z- 
component) and dephasing with the rate F2 = 7 -I- Fi/2 
(decay in xy-plane), where 7 is the pure dephasing rate. 
In the Rabi rotation, the phase ip controls the axis; e.g.. 
If — and if — Tr/2 cause the spin rotations around x 
and y axes, respectively. Combining pulses with differ- 
ent if, one can controllably rotate the spin, while \a\ < 1 
subjected to the incoherent processes. 

The atom generates two coherent waves propagating 
forward and backward with the current field 
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fiFi 



^ik\x\—iujt 



(3) 



found from the expectation value of the current operator 
P — {hVi/ (j)p)ia^ of the atom situated at a; = 0. Equa- 
tion ([3]) presents the dipole emission proportional to (cr^ ) , 



2 




(b) 



Atp 






— ► 


— ► 


► 


P 




M 




R 



r '^atom" \ 

I X XX X I 






FIG. 1: (a) Measurement circuit diagram. (b) General 
pulse sequence (upper panel) with schematic atomic dynamics 
(lower panel): The pulses P and M are used to prepare and 
manipulate the atomic states, and the emission from the atom 
is detected during the readout pulse R. 



that is, to the projection of the pseudo-spin on the xy- 
plane. Therefore, the emission amplitude is a direct mea- 
sure of superposition in the atom and reaches maximum 
(K'''^)! ^ 1/2) in a maximally superposed state. The 
real and imaginary parts of the emitted wave (0- and 90- 
degree phase in respect to the driving wave) give (ay) and 
{(Jx), respectively. In the case of the atom driven by the 
microwave with ip — from its ground state {0,0, —1}, 
the emitted wave is purely real because (ax) — 0. 

In Fig. [if a), a schematic diagram of the measurement 
setup is presented. We prepare and manipulate the atom 
by one or two sequential microwave pulses (denoted by 
P and M) of lengths Atp and Aijvi applied at times tp 
and [Fig- IHb)]- A continuous microwave is split into 
two channels and chopped independently by the rectan- 
gular pulses P and M with the choppers consisting of 
high-frequency mixers. The microwave in the P channel 
acquire a phase shift ip (=0 unless it is specified) pro- 
duced by a computer-controlled mechanical phase shifter. 
The two channels are combined into one, and the signal 
is delivered to the sample in a dilution refrigerator (at 
T w 50 mK) through a coaxial cable. 

The output signal is amplified by a cryogenic amplifier 
and a room-temperature amplifier and then chopped by 
the readout pulse R of a length Atp — 50 ns (Atpi ^ r2"^) 
at time tp after the manipulation is completed. The sig- 
nal is detected either by a vector network analyzer (VNA) 
in the homodyne measurement of the coherent emission 
or by a scalar spectrum analyzer (SA) in the power mea- 
surement. As it follows from Eqs. ((Ij and ([21), during 
the readout pulse R the current amplitude of the co- 
herent emission is /+(t) = (ftri/0p)(cr-(<R))e"^2(*"*'^) 
and the forward (backward) emitted power is P(t) = 
(;ia;ari/4)(l -I- {pz(tp)))e-^^^^-^^'^ . The signal is inte- 
grated over Atp and the pulse sequence is repeated with 
a period of = 250 ns, resulting in the averaged detected 
signals by VNA and SA /+ = hV i / ((jipV 2Tr) {ij- {tp)) and 
P hjJa/{'iTr){l + {(Jz(tp))), respectively. We will omit 
tp in further notations. 

Coherent and incoherent dynamics of the system is 
studied by measuring the emission from the atom. A 
single microwave pulse P of varied length Atp is applied 
to the atom with the readout pulse R following right after 
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FIG. 2: (a) Pulse sequence for the preparation and readout the 
pseudo-spin and the corresponding rotation in yz-plane of the 
Bloch sphere, (b) Power spectrum of the free-induction decay 
emission peak with a long excitation pulse Atp — 100 ns. (c) 
Evolution of the atomic population measured by the emitted 
power (upper panel), and the normalized dipole moment (o"~) 
derived from the coherent dipole emission (lower panel). The 
dashed lines mark the points where the emission is completely 
incoherent ((f") = 0), corresponding to the extrema of the 
emitted power (oc 1 -f (o-^)). 



P as shown in Fig. [2][a) . The SA spectrum of the inco- 
herent emission from the atom excited by a long pulse 
with Atp = 100 ns (to saturate its population to 50%) 
is shown in Fig. ^h) by the red dots. By fitting it with 
a Lorentzian (black curve), we obtain the dephasing rate 
r2/27r = 9.4 MHz. The peak power contains informa- 
tion about the atomic population. The upper panel in 
Fig. ^c) shows (iTz), measured by the emitted power in 
5 MHz bandwidth at Wq. On the other hand, the lower 
panel shows the dynamics of (cr") derived from the co- 
herent dipole emission measured by VNA. The behavior 
is consistent with the rotation of the spin in yz-plane: 
Re((icr^)) = (o'j^)/2 (blue curve) exhibits decaying os- 
cillations, while the imaginary part lm{{ia^)) = {ax)/2 
(black curve) remains nearly zero. The total emission 
becomes completely incoherent in the extrema marked 
by dashed hues, because at these points (cr~) = 0. The 
curves of Fig. [2fc) have also been used for calibration. 

The oscillations decay, while the emission satu- 
rates to a finite stationary level defined by («cr^) = 
(rifi/2)/(rir2 -l- n^) {= {cry)/2). The stationary emis- 
sion at different driving amplitudes is exemplified in 
Fig. [2Ia), where the evolutions of (cry)/2 are shown for 
n/2Tr = 140 MHz (red curve), 44 MHz (blue), and 
14 MHz (black). Here, the saturation (stationary) level 
increases with decreasing fi. As it follows from Eq. ([T]), 
the oscillations decay with the rate ri/2 -I- r2/2, because 
the rotating pseudo-spin is subjected equally to relax- 
ation in z-axis as well as dephasing in y-axis. The dashed 
black curve tracing the oscillation maxima is the expo- 
nential decay with (ri/2 + T2/2)/2t: w 13.5 MHz. 

Having full control of our system, Fi and r2 can be 
measured separately. Figures. [Sfb) and (c) present relax- 
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ation processes with decay rates T2 and Fi . To determine 
we apply a 7r/2-P-pulse [the lengths of 7r/2- and tt- 
pulses are defined from the first maximum and first zero 
on the blue curve of Fig. HJc)] which brings the pseudo- 
spin to y-axis, and then measure the emission after a 
delay AtpR. The decay of (ctj,) as a function of AtpR is 
shown in Fig. [2Ib) . To measure the energy relaxation, we 
apply TT-F-pulse for preparing the pseudo-spin up (along 
z-axis) , and after a delay AtpM , during which the popula- 
tion (1 -|- {az))/2 decays, we apply 7r/2-M-pulse rotating 
the spin from z- to y-direction, as shown in Fig. |3jc). 
By fitting the experimentally measured emissions with 
exponential decays, we obtain Fi/27r = 18.3 MHz and 
F2/27r = 9.1 MHz (F2 « Fi/2), which indicates that pure 
dephasing is negligible in the system. These numbers are 
consistent with the ones separately determined above. 

With the demonstrated state manipulations, one could 
perform in principle the quantum state and process to- 
mography, obtaining full information about the time evo- 
lution of the atom. Here, instead, we demonstrate deriva- 
tion of incoherent properties of our system by measuring 
only the coherent emission. Namely, we reconstruct the 
two-time correlation function of fluctuations, derive in- 
coherent spectrum of the resonance fluorescence (Mol- 
low) triplet and compare it with the spectrum directly 
measured by SA. To derive the correlation function the 
full set of tomography measurements is not needed. We 
propose a convenient procedure for the correlation func- 
tion derivation from experimental curves of the coherent 
dipole emission. 

The incoherent emission power spectral den- 
sity in either direction is given by >S'(a;) = 
(1/2^)Z /_!1(A/- (0) A/+ it)),,e'^'dt, where 
Ao = o — (o)ss denotes deviation (fluctuations) of 
the operator o from its steady state (o)ss = linit^oo(o(t)) 
[l9| and Z {— 50 51) is the characteristic impedance of 
the transmission line. Substituting from Eq. ^ and 
using the atomic relaxation rate Fi — hjjacj)^/ {h^ Z) into 
the two directions along the line, we arrive at 



(a) 



00 

S{u:) = / (Aa+(0)Aa-(0)sse^"*di. 



(4) 



According to the quantum regression formula [l9| 
the steady state of the two-time correlation function 
(Ao'+(0)A(T^ (t))ss can be found by solving the equation 



ds{t) 
dt 



Bs(t) 



(5) 



for s{t) = (A(T+(0)A(T(i))ss (where s = {sx,Sy,Sz}) 
with the initial conditions s(0) obtained from the steady- 
state CTss and then by calculating (ACT"'"(0)Acr~(i))ss = 



j{t)]/2. For the strong drive {Q. ^ T2) where 



x,y,z /ss 



)ss 0, the initial conditions are simplified to 



1 



1 



(6) 
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FIG. 3: (a) Evolution of (ay) /2 at diflerent driving ampli- 
tudes n/27r = 140 MHz (red curve), 44 MHz (blue curve) and 
14 MHz (black curve). The stationary value of (cr~) strongly 
depends on the driving power. The dashed black curve de- 
picts decay of the fast oscillations with rate (ri/2-f r2/2)/27r 
= 13.5 MHz. (b) Direct measurements of dephasing. The 
red dots show (ay) as a function of the delay Atpn between 
the manipulation 7r/2 and readout pulses. The black curve is 
an exponential-decay fit with r2/27r = 9.1 MHz (c) Measure- 
ments of relaxation. The measured decay of (ay) (projected 
{(Jz)) is presented by the red dots. It changes the sign when 
((Jz) becomes negative. The black curve is an exponential fit 
with ri/27r = 18.3 MHz. 



The solution for s(t) can be found from Aa(t) = a(t)—ass 
by measuring the coherent dynamics of ^(i), governed by 
Eq. ([T]) for two initial conditions a'{0) = {1,0,0} and 
a"{0) = {0, 1, 0} and subtracting CTss « ^{t) at t > F2. 

Here, however, we demonstrate an equivalent but 
slightly different approach, which is more practical ex- 
perimentally. We measure evolution of a{t) with two 
pairs of opposite initial conditions: <?(-|_)(0) = {±1,0,0} 
and «?('j-)(0) = {0, ±1,0} prepared by applying 7r/2-P- 
pulses with ip = ir ^ 7r/2 and (p = tt/2 =f 7r/2, re- 
spectively. Differences in the pairs give solutions of 
Eq. ©: s'{t) = [a[^^{t) - a[_^it)]/2 for ."'(O) = {1,0,0} 
and s"it) = " '?(-)W]/2 for s"(0) = {0,1,0}. 

The initial conditions of Eq. ^ can be rewritten as 
s(0) = [s'(O) + is"(0)]/2 and, therefore, the desired corre- 
lation function is {Aa+{0)Aa-{t))ss = {[4(0 + 
^[■^"(O ~*Sy Such a procedure provides robustness 

against small errors in ip and does not require separate 
measurements of <7ss- 

The evolutions of the x, y-components of Act' and Act" 
for ct(0) = {1,0,0} and {0,1,0} are shown in Figs. HJa) 
and (b). The real and imaginary parts of the correla- 
tion function (ACT+(0)ACT~(i))ss derived from those plots 
are shown in Fig. |3{c) by the red and blue curves. The 
dashed black curve shows the correlation function alter- 
natively calculated from Eq. ([S]) using B with Fi and F2. 
The correlation function of a superconducting quantum 
circuit has also been recently measured more directly in 
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FIG. 4: (a) and (b) Evolutions of {Aax,y{t)) with initial condi- 
tions S{0) = {1,0,0} and {0, 1,0}, respectively, (c) Two-time 
correlation function of fluctuations (A(T"'"(0)A(T~(t))as calcu- 
lated based on evolutions in (a) and (b). The red and blue are 
the real and imaginary parts, and the dashed black curve is a 
real part of the solution obtained from Eq. (5) using B with 
experimentally deduced parameters. In the experiment, the 
points at times shorter than 0.8 ns could not be measured. In- 
set: Spectrum [S{f) — 2ttS{lj)] of the resonance fluorescence 
triplet as a function of detuning Sf = {uu — LJa)/2n derived 
from the two-time correlation function according to Eq. ^ 
(black curve) and measured directly in the frequency domain 
SA (red open circles). 



Ref. [20|. By taking Fourier transformation of the cor- 
relation function [see Eq. (|4])], we obtain the spectrum 
of resonance fluorescence, which coincides well with the 
experimentally measured one as shown in the inset of 
Fig. He). 

In conclusion, we have studied time-evolution of co- 
herent and incoherent dipole emission from an artificial 
two-level atom strongly coupled to an open transmission 
line. The correspondence between coherent and incoher- 
ent emission has been shown. We demonstrate derivation 
of the two-time correlation function and reconstruction of 
the spectrum of resonance fluorescence (MoUow) triplet 
from time-dependence of the coherent emission. 
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